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A study of the internal structure of the neutral she&t in 
t,he g c oina gn e t ic .tail has been made from data obtained by the 
. NASA-GSF C magnetic field experiment on the Explorer 34 space- 
crart during its tail passage in the. first half of 1968. The 
data,. used in the analysis are individual measurements of 'the 1 
vector magnetic, field at 2.56 see. intervals. The experimental 
results consist or statistical studies of field orientation and th< 
component as a function of field magnitude. The results 
do not support nearly one-dimensional field models with char- 

t , , • * * * l 

acw eristic lengths for fi eld variation parallel to the neutral’’ 

sheet much larger than the neutral sheet width. The principal 4 

conclusion from tne data points toward consistency with a cuasi- 
> ' - 

-periodic (possibly turbulent) structure with a tendency to 
r-ormation of magnetic loops as one might expect from stability • . 
studies. Numerical results, obtained, using a simple model with 
a periodic field pattern’ in the field reversal region are com- ' ‘ 
pared with the observed internal structure. ^ ■ ' 
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7. ho presence o£ ( a ’ neutral sheet in the tail of the gcoinagnoto- 
spherc wae first established by Ness (1965). Later, other 
measurement s confirmed the neutral sheet as a permanent 
large-scale structural element of the goomagnotosphore (e.g. 
review by Ness, 1969). For instance, Bohannon (1968) has 
shown that the neutral sheet extends well beyond the orbit of 
the moon. The picture of a plane current sheet separating 
regions of magnetic field with opposite directions has become 
useful for the construction of quantitative models for the 

, i 

gross ‘ma gnetospher ic ficild. (See review by Roederer (1969 ) . ) 
Sparser and Ness (1967) concluded from IMP-1 data taken at 
radial distances between 15 -and 32 R.~, that there is an average 
south-north field component of a strength of 1-4 gammas super- 
imposed.; this points towards some connection of the pulled- ‘ 

-out geomagnetic field lines through the neutral sheet. In . 
this picture the field is' still approximately one-dimensional, ‘ 
in that the s-dependence dominates, 

Such one-dimensional field descriptions are valid if one is 
interested in the average .field structure on time scales 
larger than several minutes. On t smaller time scales one 
fir.-ds a more complicated structure, to' which we shall refe’r' 
as the "internal structure of the neutral sheet 1 '. '.** . 

•i * . * 

Evidence for ‘the existence of the internal neutral sheet • • • 

. * 

structure was provided by Mihalov et al. (1968 ) and by Mihalov,. 
et al, (1970), who discussed the occurrence of both nositive 

i • * 

and negative values of B (solar magneto spheric coordinates) 
in neutral sheet traversals of Explorers 33 and 35. In -ad- 

* ■ ■ ■ g 

ditior. to the picture of geomagnetic field lines pulled out 

. I. ■. ... 

into the tail, the formation of closed loons was mentioned as 

. : .• • . - ■ ; ; ■ jf’ » • 

a possibility, \ /. f...---;;. : f ft. if pi/ '‘//t: f 

In this note we investigate the internal neutral sheet struc- ./,/.. 
ture op the basis of magnetometer data from Explorer 34, using the 
best time .resolution available , It is our aim to collect more 
quantitative information and to discuss the bearing our findings 
have on existing theoretical models . The analysis leads to a 
tentative concept of the topology of the structure.’ 
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The NASA-O.SFC magnetometer on the Explorer 34 spacecraft has 
bean described elsewhere (Fairfield* 1969 ) . 

here we only state that the magnetic field vector is measured 
at 2.56 second intervals with a precision 1 of ± 0.16 y and ’an 
accuracy of 0.1, y; the equivalent RMS noise in the 0-5 Hz pass 
band of the magnetometer is 0.028 y. 


In February-March '1968- the orbit of Explorer 34 was in a- po- 
sition favourable fo.r ,,t‘he present purpose. Neutral sheet cros- 

i 

sings were observed on orbits 61 to 66 in the region 


27 R_ % 
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13 time intervals have been selected, giving a total of 16351 
samples of the magnetic field vector (Table 1). During these 
intervals there was an enhanced tendency for F = -|bJ to assume 
relatively small 'values . Figure 1 shows the corresponding orbi 


sections* 


•The average z position of the satellite during each of the 

u ■‘ilR 4 ” • 


13’ intervals studied is included in Table '1. It is seen, in 

* * , # 

the table .and Figure l„that the observed neutral sheet position 
most frequently, was within 1 R £ of the solar magn etospher ic 
eauatorial plane (z = 0)'. Three of the intervals selected 
occurred when* the neutral sheet position' was -5.3, -4.6 and ■ 
*3.6 R„ from 'the equatorial plane. Such departures are not 
unexpected at this radial distance from the earth, ^32 R^ , 
when account is taken of the tilting of the • geomagnetic dipole 
utilizing formulae such as given by Russell and Brody ,(1967). 


The typical behaviour of the magnet ic 1 field during the' neutral 
sheet traversal is shown in Figure 2. In Part a, each data 
point represents the time average .over 20.43 seconds of eight 
individual measurements obtained at 2.56, second interval^’. Th.e 
plot of data -in Figure 2a is 'typical of that employed for . : • 

selecting the more restricted intervals of neutral sheet cros- 
sings used in this paper. Tho time interval selected in this 
particular data set is indicated by the, hatched markings.’ 


s 


In Part b arc presented the detailed plots of data obtained 

at 3. 50 second intervals which typify the data sets utilised in 

this analysis. It is clear why high -resolution data arc ncc- 
» 

ossary in studying the detailed structure of neutral, sheet re- 
* 

vorsal regions." With very low sampling rates or through use 
of averages , such as employed in the very early exploratory 
studies in the IMP-1 work, much detail is lost. 


Particularly noteworthy in those figures is the very rapid 
change of the magnetic field observed at the position of 


the satellite. Two possible explanations exist. One is that 
the spacecraft is sampling a .time -stationary but spatially 
variable structure of tfye neutral sheet ; the other that the 
fixed spatial structure of the neutral sheet is moving back • 
and forth relative to the satellite position. There is no 
unique way in which these single satellite data can be studied 
to determine which of the two possibilities dominates. 


It is in this spirit. that a statistical study of data ob- , 
tained in the neutral sheet reversal region was considered 
» as one • possible way to determine which, if any, of the pos- 
sible models for the neutral sheet structure were more, con- 

• , • * , 

• si stent with the, observed data. The Idea was to compare 

* * * * 
the statistics of a time- stationary spatially varying theo- 

• i 

retical .structure with the observed . data , assuming that the, 

* i * . , 

essentially spatially 'fixed satellite was then sampling •( in 
time) the spatial structure of the field reversal region. 


The number oi: sample points utilized in each of the 13 data 
subsets varies considerably (see Table 1)’. A consideration of the 
state of the magnetosphere is also relevant for possible, 
non-typical conditions. The planetary magnetic activity index 
Kp for February 22 was 3- 3* 3 3 2-3- 1 0+ while February 25, . 
was, one of the five quietest days in February, Kp being 1 It 


0t 0 1 0 n* 1 2 2 1 


On February 27., one of. the five quiet - . days,, 


it was 2 t 1- 1 1- ,1 2- 4- 3. Thus for the major portion of 
•the data studied we can reasonably expect the state of the' 


~ 4 - 


magnetosphere to bo typical of t h o quiet magnetosphere , with rut 
significant distortion duo to solar activity. 


We now turn to a discussion of the data in some more detail. 
We have grouped tho material in intervals according to the 

in/b^T 


value- of E 


F 


Table 2 gives the intervals 


V „ o 

•a tB . 

x y a 

and the quantities of major interest for the present discussion 


An important aspect of the data is the fact that there is ? 

• . ■ ' ■ * 

'•ubstantial number of cases* wh ere F is small:- F < 2y in 
about 13% of the samples, 1 . 


As will bvuome apparent in the following section , is 

• » ’ +4 » 


'more interesting for* the comparison with existing model; 


than 3^ and thus we concentrate on the properties of . 


'he average value of B r , B t , is rather small for small F 


’and then increases to 'v 2y, the overall average, being - 

* * 

1.68 y» in fair agreement with earlier results (Speiser and 

* 

Ness j’ i967 ) ‘ for essentially .the same region of the tail. 


A feature which is particularly important for the present 
purpose is the fact that there is a .significant fraction . 


of cases • (24% )• with negative values of 3 . In ' the ' first 

z 


interval (0<F<1 . 25 y) about half of the field vectors are 

.pointing southward. The inclination of the field lines with 

\ ’ ' ■*' : 

the x,y-plane » measured by J 0 1 , is fairly large (27°) 
in -the .first interval, and even increases with F reaching a 

i * 

• « 
maxim'um in' the 3rd interval. . ' 


It is interesting to‘ note that these measurements are car- 


ried out at geocentric distances smaller than 35 R_ , where 

♦ v ' 


Mihalov at al . • ( 1968 ) found essentially (one exception) no 

case of negative .B .• .The reason may well be the higher time 

. ■■■•#'. , * , 

resolution • , seeing that there is still considerable 

point to point variation at the present resolution. ' 


■" ; 


- 5 - 




“s tr:: 


rtt .?••••# • • * *** - «?■!*•* -***# r- 


» 


We new approach the problem of the internal structure by 
discussing what one would expect in a simple ono-dimonsional 
picture where the field lines arc pulled-out dipolo line, o’. 

This model seems to bo consistent with observations on suf- 
ficiently large space and time scales (Fairfield and Ness » 

1970). 

* * i 

Since the present observations are confined to a rather 

* • 
small radial distance interval , wo can ignore spatial dependence 

along the tail and conclude that in such a simple situation 

B must be approximately constant (V*B=0). In this case all 

•i , 

data points would lie on a‘ horizontal line in the B , B 

z * x 

plane (Figu’re 5). This model is .however not adequate because 

i • • * * , 

it does not take into account the fact that the neutral sheet 
moves such that the instantaneous direction of the sheet nor- 
mal tilts away from its average orientation. (Mihalov et al . , 

1970). , ' 

» 

If we interpret this tilt in terms, of a wave .with sufficiently 
large wavelength we. would expect that every magnetic field 
vector at a given point in space oscillates around- its average 
orientation with 1 a fixed angular amplitude. In ‘that case the 
data points would lie in a wedge-shaped region such as shown in 
'figure 3. (The boundaries correspond to horizontals in suitably 
rotated coordinate systems, ) 


Figure 4,. part a, shows the observed field orientation, 
in 'each magnitude interval separately. Since we are dealing 
with the projection of thejfield, onto the * sus-ulane we -have 

, * : /r 1 1 * /— ~ 

now ordered the measurements with respect to G = >/ 3 t> 
using the same intervals as for F (Table 2). 


~ I 2 
+ B 
X z 


There are two distinctly different regions. For .0 X £ Y the- • 
picture is roughly consistent with * that given in Figure 3 . .> 
From ,fhe established- gross picture of the neutral she.et it 
seems reasonable to visualize regions of larger G as lying 
on the average, further away from the center of the neutral 
sheet than those where G is .smaller . We therefore conclude 


6 


*t**r ' *ri~ z ■ •? 


•* r*"r- 


=5 *■ 


that in the region whore G £ 5y, the local tail .struc- 
ture can be explained by a one-dimensional field configuration 
with “i rt £‘2.3y and with long- scale waves superimposed which 

M * " , 

lead to an angular distribution of half -width ^20° (straight 
lines in figure 4 a drawn by hand have a elope of 22°). In 
terms of A If van wavus this corresponds to an amplitude of 
* 7y in the full tail field of 20y. This orientatio’n effect 
then' would lead to the B < 0 observations for G £ 5y. 

•St 

v * 

Xn the regime G < By this explanation clearly fails to bo 
applicable. Instead of being constant B drops to aero 

■ » » 52 * 

* S *" - r . 

for decreasing values of G, and instead of > 0 in all 

.eases, the first G-interval contains about 49% of data points 
with B < 0. Table 3 gives some, of the relevant properties 
for G £ 5y. Obviously the field is considerably weaker than 
expected from a one-dimensional picture and is more ran- 

•r . * 

domly distributed. 1 • In the first interval the distinction 
between north’ and south has almost completely disappeared. 

* 

3. Interpretation 

• t k 

Several models exist which depart from the concept of a one- 

- * f * » 

-dimensional tail field. ‘The reconnection model has as an 
essential feature the presence of a neutral point somewhere 
in the tail (Dungey , 1961 ; Axford et al. , 1965): 


Ppssler and Hill (1970) have suggested a model which contains 

a neutral point -relatively near the earth, e.g. at.lSR^ radial 

■ distance, and a neutral sheet with a weak negative B -component 

2* 

beyond. 


A quite different approach to a two-dimensional structure of 
the tail is relat.ed to a stability analysis of a one-dimen- 
sional neutral sheet (Furth , 1962 j . Pfirsch , 1962 ). . •» . 


-a 




For a, large variety of circumstances (Schindler, lSGGi 
Schindler and Soop, I960) one finds instability against 
tearing, i.a. the magnet ic field linos break up and form 
closed loops. The growth rate was computed by Laval ot al. 
(1966). Those results ‘iwere first applied to the geomagnetic 
tail by Coppi et al . (1966). Bis’kamp ei al. (1970) and 
Biskamp and Schindler (1971) have recently studied problems 
related to the non-linear evolution of this instability. 
Linear theory predicts stabilization by anisotropy (Laval , 
and Peilat, 1968), It is however not difficult to visua- 

» 

lize non-linear perturbations which are able to grow. The 
tearing mode might also be driven by the presence of a 
(collective ) * resistivity (Furth et al., 1963 ). The question 
how a magnetic field component normal to the neutral sheet 
influences linear and non-linear stability is essentially , 
open. 


If one takes the reconnection model literally, one expects 
one neutral point. < The only possibility to reconcile this 

i ' " . * 

picture with the present results is to assume that the 

, « 

neutral point position rapidly oscillates. 

* • t 

i 

In the case of the present Explorer 34 data, the neutral 
sheet observations occur when the satellite is in the dis- 
tance range of 30 to 34 earth radii. Thus a rather local' struc 
ture is being sampled' by this spacecraft. The data published ' 

•V i . -• 

by Kihalov et al. (1968) alee contain both '.northward and "• 

, 1 i • • ’ $ 

southward pointing vectors, for instance at radial distances 
exceeding 70 R„. This would suggest an amplitude of at least 

t . * . . *** * • « 

40 Rp . . ’ ’ 


In order to obtain a period which one would have to associate 
with the oscillations of a single neutral point, we have com- ’ 

puted the average 'time between sign reversals o"f B . This 

V *' 2 ' Z 

time is of the order of 10 
* 

shown in figure 2b contains 


sec . 

For 

instance , * 

4S B 

i 

,-sign 

reversals 

art • 

The 

same order 

ould 

folio 

M from the 


' ' \\ * ■ • 

by Kihalov et al. (1-OS8), who found a considerable reduce io: 


V 5 


in uef,r*\ivQ-?„ observations if they looked far* ca nr , 3 whoV'j 

•a * 

was email (F < 2) over about 100 sue., averaging over* that 
time interval . For* a sinusoidal neutral point motion this 

would *ivo a lower bound or. the velocity amplitude of about 

9 >• 

10 cm/sec, Since we arc using data collected on three days 

only, we cannot safely say whether t the structure observed is 
typical for all times* On the other hand, even if our cofi- 
mate of 40 R.„ for the amplitude is too high by a factor of 

ti 

10, wo still would find a velocity which exceeds 5 times the 

Alfv&n velocity. There .does not seem to exist evidence for 

* * » 

such a supersonic oscillation. 


‘Furthermore, it is very difficult, if not impossible', to re- 
concile the picture of a single oscillating neutral point 
with the fact that for *>/b * £ 5y there seems to exist 

a simple one-dimensional structure with an average positive 
normal component, which at radial distances of about 00 R„ 
is of the, order of 2 y . 


.The, second type of structure discussed above, i.e, local 

» * * ■* 

concentrations of the electric current, in -the simplest ; 
picture giving rise to magnetic loops is not. subject to 

t * ’ 

these difficulties. Since it allows for many neutral points, 

’ 

the rapid change of sign of B is easily explained by a 

convective motion of the current concentrations . 

*. - ■ . . • ■■ * 

< ■ , > . * - • 

Clearly, the presence of current concentrations, whose ’char- 

* 

act eristic dimensions are expected to be roughly of the . order* 

of the neutral sheet width, is not sufficient to explain all 

, * 

features of the observations. One has to add a net magnetic 
flux in the south-north direction. ■ 

. 1 - . . . if ■ . f 

A simple model combining the two features is discussed in 
the following -section . ‘ 

, ■ • • » . . * 

4‘. A simple' model , r 

Simple .models which have the desired topo logy in the x-z 
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’H* 


plane can bu conveniently represented by a suitable fern 

of the y -component of the vector potential, A(x,n). A 

possible choice is 

. • 

K 1 


A(x,») = B L 
o 


cos 


cosn 


K I 

$ 


■log(cosh f t a ; &*;) t b p ( 

'■ JiS h 

L J 


where B q is the field magnitude away from the shoot and L 
,is its width,, a is a measure for the size of the loops 
and b3 o is the magnitude of the superimposed homogeneous 
field in the ^"direction, For D - b = 0, a « 1 this model 
reduces to the linear tearing mode in the limit k L->1 , where 

* JK 

. k is the wave number in the x-direeticn. Choosing D i 0 

provides the possibility of confining the mode to a smaller 

> 

region around x = 0 , • 

This model is not self-consistent because we cannot explicit- 
* * 

ly write down particle distribution functions which are 

in equilibrium , with the fields chosen (this is only possible 

for *b, D s 0J. For the present purpose, where we are mostly 

concerned with, the qualitative structure of the field, this 

is however not a serious drawback . M. Soop (private communis-' 

ation) numerically obtained field patterns which solve the 

steady state equations approximately. 1 ' In some of the cases 
# 

studied the magnetic field has the same topology as the simple 
model given above. 

A possibly more serious deficiency of the present approach 

* - » 

is the omission of the B -component. We will return to this 

■ , . y 

point in the following section. 

• ■ ... • . • 

* „ . t 

Choosing B q = 20y and b = 0.1 (Behannon, 1968) we have es- 
sentially two parameters to, vary, a and D. It turns out 
that the case a - O.S> and ,D - 3 is able to reproduce the 
essential features of the observations reasonably well, ; Be- 
cause of the arbitrariness of our model we do not, aim for 
excellent quantitative agreement. There is probably a different 
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Such 


iiei* of values a, D which give on even bettor fit, 
an optimisation was not done because good quantitative 
agreement in not essential for the conclusions wo can draw 

from our ad hoc model. 

« * * 

‘ * » 

Figure 5 shows the field lines for the case choson. Since v/c* 
have' not added ne'utral sheet tilting in our model and since 
their effect does not seom to be essential for B J 5y (figure 
4) we have ignored orientation effects for 6 ^ 5y and drawn a 
wedge-shaped region for G X 5 Y • T,here seems to be good quali- 
tative agreement* A comparison between parts a and b of 
figure 4 shows that the general behaviour is reproduced. 


* 

So far we have discussed only moments of the probability 
distributions ‘involved , It is also interesting to compare 
the distributions of the magnetic field components themselves, 
..Figure -6 gives the observational and model distributions of B 

• < , i 54 

for* the first 6 intervals. Although the corresponding' distri* 
buttons show considerable deviations they have a number of 

i * 

qualitative properties in common. , The main features are the 

, * 

.occurrence of two maxima, the larger peak corresponding to 

* . * 

a larger value of . It seems not very surprising that the 
experimental curves are considerably smoother than those of 
the model. This may be mainly due tc the presence of .a more 

irregular structure than the one discussed here, to super- ' ’ 

1 

imposed large-scale waves, or to averaging over a number of 

• f i * 

different configurations. * * . 


. \ . . . ■, : . # * ♦ : 

In this analysis we have concentrated on the pro j ection '*of B 

into t the x , z-plahe . 'Although the behaviour of seems to 
be less important for the present purpose, it is interesting 

to note some qualitative properties. Comparison between Tables 

* ■ ... 

• 2 and 3 shows no marked changes when going from F to G . as 
the ordering quantity. ‘ 


Furthermore , the distributions 
pact to F ) look ^similar . This 


of B and B (ordered with res- 

■ y 3 

suggests that the internal struc 


E*S S * 


Is three -dimensional . Such field configurations are qualita- 

* 

tivcly consistent with non-linear single mode growth of the 

i ' * 

tearing instabilities where one expects the pinch instability 

to become important (Biskamp et al. (1970)). 

» 

5. Discussion ’ • 

♦ » 

* ♦ 

* 

We have seen that for the average tail configuration asso- 
ciated with the* time intervals during which the observat lens 
discussed were made , a simple two-dimensional field line model 


can be apul’ied for 


+ B ' > 5y. l’n the regions 
* z 


✓ 2 2 

3 +3 < 5y there is considerable deviation from this 

x a , 

picture. The essential feature is the occurrence , of region 
with very weak magnetic field which we interpret as’ neutral 
points. The concept of one neutral point which oscillates 
along the tail is not supported by the present observations. 

Therefore we have investigated the consequences of a model 

* . t 

with many neutral points. A rather simple model of that 

• , . * , 

type has a number of qualitative properties in common with 
the observational material. In particular, the occurrence 

* f i , . * 

of negative values of , and the average and spread of the' 
field orientation, are quite -well reproduced (Fig. 4).. We 
therefore conclude that the field topology given in figure 
5 is consistent with the observations discussed. We do riot, 
of course, suggest any detailed agreement nor do we claim 
that the current concentrations are stationary., In fact, one 
expects that -they are convected along t.he tail, and that they 
are considerably more irregular,. As a- possible source of the 
current concentrations we have discussed 'the tearing instab- 
ility, which might be driven either by collision-free proces- 
ses or by a collective resistivity. . 

An- alternative process possibly leading to the change "of field 
structure at *v5y -is" that magnetosheath plasma may have access 
to the neutral .sheet in a relatively thin layer; . 


‘ - 12 - 


* t * 


disturbed magnotoshoat h finin n < 

b ncatn l3>cld l3 * nc s may connect onto the 
tail field in a random wav rivancr *>■}« 

This „ 7 3 E riso to ,na "y noutral points. 

hlS , PP0 ° eSS " Xsht lead *° a ^ B netic topology similar 
to that produced by turbulent tearing. 

♦ , 

Models irf which the magnetic field' close to the neutral 

1S £, ™ 001;l1 ^ 1,e> vat, ying on scale lengths of the tail 
dimeter or tail length) are faced with the difficulty that 

r a *J L * V xere , L ls the distance over which the quantity 


' r !* I ^ rea0hes ^ (Fi S- “a). If for instance, we choose 
L 1000 km, we find from V-B = 0 that ^2- + Hx ^ 1°Y . . . 

cannot be reconciled with smooth field ^riatlV ^ 

■ ’ f 

, I , 

configurations containing, instead of a series of. neutral 

points, a neutral surface which may be disturbed by waves 
woula not predict the net magnetic flu* iri z-direction ob- 
served on time scales long compared with the wave, period, 

. The- fact that one-neutral-point configurations seem, unlikely 
h-v e t h PP X ° able ' doeS not mean that magnetosphere models, which 
'do ions" ‘I?™'” their main pre- 

dels th‘ ^ tUrn ° Ut th3t at lEaSt in some « the mo- 

the gross properties of the magnetosphere • are still the' • 
same. A detailed reconsideration seems however, necessary. ' 
In- particular , -for processes on length scales of the' order 
of the neutral sheet width and on time- scales less than 'a" - 

ocructures seem to play an important role, 
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1 *: 


Figure 2 ; Typical neutral sheet crossing obsorved by the mag- 
netic field experiment on Explorer 34, In Part a, 
the solar magnctospheric latitude (0), longitude 
(0),*and field magnitude (F) are shown for sequence- 
-averaged data (see text).' A sub-interval of the 
data on an expanded time scale and finer time re- 
solution is shown in part b (see text). 


figure 


One-dimensional model with long-wavelength per- 
turbations i data points .lie approximately in wedge- 
•shaped region. 


Figure 4 : ' a ) Observed field orientations. Heavy lines correspond 
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b) Field orientations using model field (see text). 


Figure 5: Model field lines. 
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Figure 6 : Distribution functions of B ( 


a) observed > 


b ) model. 
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